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ABSTRACT 

Aims. Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. 
Methods. We observed the massive (-800 M©) starless gas clump IRDC 18310-4 with the Plateau de Bure Interferometer (PdBI) at 
sub-arcsecond resolution in the 1.07 mm continuum and N2H + (3-2) line emission. 

Results. Zooming from a single-dish low-resolution map to previous 3 mm PdBI data, and now the new 1.07 mm continuum observa¬ 
tions, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be 
roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the 
typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous 
initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70pm, we iden¬ 
tify weak 70pm emission toward one core with a comparably low luminosity of -16L©, re-enforcing the general youth of the region. 
The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the indi¬ 
vidual components (in the 0.3 to 1.0 km s -1 regime). Based on single-dish C ls O(2-l) data we estimate a low virial-to-gas-mass ratio 
< 0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results 
in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear 
to have very low levels of internal turbulence. 

Key words. Stars: formation - Stars: early-type - Techniques: spectroscopic - Stars: individual: IRDC 18310-4 - ISM: clouds - ISM: 
kinematics and dynamics 


1. Introduction 

Independent of the various formation scenarios for high-mass 
stars that are discussed extensively in the literature (e.g ., 
IZinnecker & Yorkel120071: iBeuther et al.ll2007l: iTan et al.ll20l4) . 

the initial conditions required to allow high-mass star forma¬ 
tion at all are still poorly characterized. The initial debate 
even ranged around the question whether high-mass starless 
gas clumps should exist at all, or whether the collapse of mas¬ 
sive gas clumps starts immediately witho ut any clear starless 
phase in the high-mass regime (e.g. jMotte et al.112007) . Recent 
studies indicate that the time span during which massive gas 
clumps exist without embedded star formation is relatively short 
(on the order of 50000 yrs ), but nevertheless, high-mass starless 


gas clumps do exist (e.g., 

Russeil et al. 201(J Tackenberg et al. 

2012; Csengeri et al. 2014 

). Current questions in that field are: 


are high-mass gas clumps dominated by a single fragment 
or do we witness s trong fragmentation during earliest evolu - 
tionary stages (e.g.. iBontemps et alJl2010t Izhang et aT] l2015h? 
What are the kinematic properties of the gas dDobbs et all 


* Based on observations carried out with the IRAM Plateau de Bure 
Interferometer. IRAM is supported by INSU/CNRS (France), MPG 
(Germany) and IGN (Spain). The data are available in electronic form 
at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 


2014) ? Are the clumps sub- or super-virial ( Tan et al. . 2014 )? 
Do we see streaming motions indicative of turbulent flows (e.g., 


Bergin et al. 2004; Vazquez-Semadeni et al. 2006; Heitsch et al 



Baneriee et al] 120091: iMotte et al.l 120141: iDobbs e~ 


To address such questions, one needs to investigate the earli¬ 
est evolutionary stages prior to the existence of embedded heat¬ 
ing and outflow sources that could quickly destroy any sig¬ 
natures of the early kinematic and fragmentation properties. 
Furthermore, high-spatial resolution is mandatory to resolve the 
important sub-structures at typical distances of several kpc. To 
get a taste of the important scales, we can estimate typical Jeans 
lengths for high-mass star-forming regions: for example, average 
densities p of high-mass star-forming gas clumps on ~ 0.5 pc 
scale in the 10 5 cm -3 density regime at low temperatures of 
15 K result in typical Jeans fragmentation scales of -10000 AU. 
Going to smaller spatial scales (< 0.1 pc), the embedded cores 
have higher densities in the 10 6 cm -3 regime that result in much 
smaller J eans fragmentation scales on the order of 4000 AU (e.g., 
iBeuther et al.ll2013bb . 


Our target region IRDC 18310-4 is a 70yum dark high-mass 
starless region at 4.9 kpc distance that was previously observed 
with the PdBI in the CD configuration in the 3 mm contin¬ 
uum and the N 2 H + (l- 0 ) emission with a spatial resolution of 
4.3" x 3.0" (iBeuther et all 2013 b). We identified fragmenta- 
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tion and core formation of the still starless gas clump, and ap¬ 
proximate separations between the separate cores are on the or¬ 
der of 25000 AU. Fragmentation of these cores could not be 
resolved by the previous observations. The simultaneously ob¬ 
served N 2 H + (l- 0 ) data exhibit at least two velocity components. 
Such multi-velocity-components are the expected signatures of 
collapsing and fragmenting gas-clumps (ISmith et all12013 ). 

In a hierarchically structured star-forming region, the cores 
are the supposed entities where further fragmentation should 
take place and which likely contain bound multiple systems. 
Studying this region now at 1.07 mm in the dust continuum and 
N2H + (3 - 2) at sub-arcsecond resolution reveals the hierarchical 
fragmentation and kinematic properties of the densest cores at 
the onset of high-mass star formation. 


tion of the frames pointin g product based on the Herschel gy¬ 
roscope house-keeping (ISanchez-Portal et all 2014b . This often 
improves the absolute pointing accuracy and mitigates the point¬ 
ing jitter effect on individual frames. We used these corrected 
frames and performed the individual detector pixel distortion 
correction, the de-striping, the removal of_l/f noise^and the final 
map projection using Scanamorphos ( Rousse l 2013 k version 
24, employing the “galactic” option. The resulting full width 
half maximum (FWHM) of the point spread function (PSF) at 
70, 100, and 160 jum are 5.6", 6.8" and 11.3", respectively. 


3. Results 


2. Observations 

2.1. Plateau de Bure Interferometer 

The target region IRDC 18310-4 was observed with the Plateau 
de Bure Interferometer at 1.07 mm wavelength in the B and 
C configuration in 4 tracks in March and November 2013. 
The projected baselines extended to approximately 420 m. The 
1mm receivers were tuned to 279.512 GHz in the lower side¬ 
band. At the given wavelength, the FWHM of the primary 
beam is approximately 18", and we used a small 2-field mo¬ 
saic to cover our region of interest. Phase and amplitude cali¬ 
bration was conducted with regular observations of the quasars 
2013+370 and 1749+096. The absolute flux level was calibrated 
with MWC349 and bandpass calibration was done with 3C84. 
We estimate the final flux accuracy to be correct to within ~ 
15%. The phase reference center is RA (J2000.0) 18:33:39.532 
and Dec (J2000.0) -08:21:09.60, and the velocity of rest vi sr is 
~86.5kms -1 . While the channel spacing of the correlator was 
0.084 km s -1 ,to increase the signal-to-noise ratio we imaged the 
data with 0.2 km s -1 spectral resolution. For the continuum and 
line data we used different weighting systems between natural 
and uniform weighting to improve the signal-to-noise ratio in 
particular for the line data. The resulting synthesized beam for 
the continuum and N2H + (3-2) data were 0.6" x 0.49" (with a 
position angle of 11°) and 1.13" x 0.5" (with a position angle 
of 13°), respectively. The corresponding lcr rms for the two data 
sets is ~0.6 mJy beam -1 and -15 mJy beam -1 measured in a line- 
free channel of 0.2 km s -1 for the line data. 


2.2. Herschel PACS 


The Herschel data were already presented in iBeuther et al.l 
( 12013a ). However, they were fully re-calibrated for this analysis 
with particular emphasis to applying the most recent pointing 
and astrometry solutions. The region containing IRDC 18310- 
4 wa s observed with the Herschel spacecraft (Pilbratt et all 
2010 ) within the Key Project EPoS dRagan et all 1201 2b . The 
related observations utilizing the bolometer cameras of PACS 
dPoglitsch et all 1201 Oh took place on April 19, 2011 (opera¬ 
tional day 705) and were contained in the observational IDs 
1342219060-63. The PACS prime mode was employed with a 
nominal scanning velocity of 20" second -1 . Data were taken in 
all three PACS filters (70,100, and 160 yum) and resulted in maps 
with roughly 10 arcmin field-of-view. We re-reduced the existing 
Herschel archive data, starting from the Level 1 data. These we 
retrieved from the previous bulk processing with the HIPE 12.1 
software version, currently contained in the HSA. We included 
the newly developed data reduction step calc Attitude: a correc¬ 


3.1. Millimeter continuum emission 


Figure [Jp resents a compilation of the 70 um and 1 .2 mm single¬ 
dish data dRagan et al 1120121: IBeuther et alll2002l). the previ ous 
3.2 mm PdBI continuum observations (iBeuther et alll2013bh as 
well as our new 1.07 mm PdBI continuum maps. The three 
datasets cover a broad range of resolution elements from 11" 
in the MAMBO 1.2mm data, to 4.3" x 3.0" in the 3.2 mm PdBI 
data, to 0.6" x 0.49" in our new 1.07 mm observations. At the 
given distance of 4.9 kpc, this corresponds to linear resolution 
elements of - 54000 AU, - 18000AU and - 2700 AU, re¬ 
spectively. With the different resolution elements, we can clearly 
identify hierarchical fragmentation on the different scales of our 
observations: the large-scale single-dish gas clump with a total 
mass of - 80 0 M© fragments into 4 cores with masses between 
18 and 36 M© (iBeuther et alll2013bh . and these cores again frag¬ 
ment into smaller sub-structures in our new 1.07 mm contin¬ 
uum data. For our following analysis, we only consider 1.07 mm 
sources that are detected at a > 4 cr level. We can estimate 
the masses and column densities of thes^smallest-scale struc¬ 
tures with similar assumptions as taken in lBeuther et alJ ( 2013b ), 
namely optically thin dust emission at a low temperature of 15 K 
(see also sect. 13.2b , with a gas-to-dust ratio of 150 (Draine, 
l201l[) and dust properties discussed in iOssenkopf & Henning! 
( 1994 ) for thin ice mantles at densities of 10 5 cm 3 (/q. 07 mm ~ 
0.95 cm 2 g -1 ). Main uncertainties for the mass and column den¬ 
sity estimates are the applied dust model and the assumed tem¬ 
perature. Based on this, we estimate an accuracy within a factor 
2 for these parameters. This way, we estimate the masses of the 
smallest-scale substructures to values between 2.2 and 19.3 M© 
(Table [I]). The difference of the sum of the core masses com¬ 
pared to the large-scale single-dish mass of - 800 M© is mainly 
caused by the spatial filtering of the interferometer that traces 
only the densest inner cores and not the envelope anymore. The 
peak column densities are very high in the regime of 10 24 cm -2 , 
corresponding to visual extinctions of - 1000 mag. 


In addition to the masses and column densities, the data 
allow us to estimate the fragmentation properties of the sub¬ 
sources, and we find projected separations between the newly 
identified structures between 2600 and 12000AU (Table |2]). 
Since the absolute peak positions can be determined to high ac¬ 
curacy (jm_ with the resolution 0 and the signal-to-noise ratio 
S/N, e.g., IReid et al.ill988l) . and the intereferometric positional 
accuracy depends mainly on the position of the used quasar and 
its gain calibration solution (within 0.1"), the projected separa¬ 
tion is estimated to be accurate within 0.1" or 490 AU. A frag¬ 
mentation discussion is presented in section |4~T1 


2 













































H. Beuther et al.: Hierarchical fragmentation and collapse signatures 


V3 o 
c o 
o 

O CN 
• o 

O 00 
O I 

o 


c 



o 

CM 


O 

0 

Q 


CN 

o 

00 


1.07mm continuum PdBI 1.07mm continuum PdBI 



R.A. (in J2000.0 and ") 


Fig. 1. Compilation of the continuum images in IRDC 18310-4. The bottom-left and middle panels show in color the Herschel 
70/mi image with a stretch going dark for low values dBeuther et al.ll2013bl linear and log stretch for clarity, respectively). The 
white contours in the left panel show the 1.2 mm MAMBO continuum observations starting from 4 cr and continuing in l<x steps 
with a 1<t value of 13 mJy beam -1 . The contours in the bottom-middle panels presents the 3 mm continuum data from lBeuther et ali 
( 2013b ) starting from 3 cr and continuing in 2 cr steps with a 1 cr value of 0.08 mJy beam -1 . The bottom-right panel then shows 
the new 1.07 mm continuum observations starting from 3 cr and continuing in l<x steps with a l<x value of 0.6 mJy beam -1 . The 
two top panels show zooms of the 1.07 continuum data with the same contour levels but a different color stretch to highlight the 
sub-structures. Each panel presents scale-bars and the corresponding synthesized beams. The 70yum PSF is 5.6". 


Table 1. Fluxes, masses and column densities 


S int 

(mJy) 

^peak 

( ) 

V beam / 

M 

(Mo) 

/Vh 2 

(10 24 cm -2 ) 

mml 

12.0 1 


9.7 1 


mml a 

8.2 

4.1 

6.6 

1.1 

mmlb 

3.8 

2.9 

3.1 

0.8 

mm2 

17.0 2 


13.7 2 


mm2a 

4.5 

3.8 

3.6 

1.0 

mm 2b 

4.7 

3.8 

3.8 

1.0 

mm 2c 

3.3 

3.3 

2.7 

0.9 

mm3 a 

24.0 

14.0 

19.3 

3.8 

mm3b 

2.7 

2.7 

2.2 

0.7 

mm3c 

3.5 

3.5 

2.8 

0.9 


Uncertainties for S\ nt and S pPa k are ^15% (42t. Uncertainties forM and 
Nu 2 are approximately a factor 2 ( 33.11) . 

1 Integrated over mml & mm la 

2 Integrated over mm2a, mm2b & mm2c 


Table 2. Projected nearest neighbor separations 



o 

(AU) 

3mm data 



mml-mm2 

2.9 

14200 

mml-mm3 

7.1 

34900 

1mm data 



mmla-mmlb 

0.9 

4400 

mm2a-mm2b 

0.5 

2600 

mm2a-mm2c 

0.7 

3500 

mm2b-mm2c 

1.1 

5500 

mm3a-mm3b 

1.5 

7600 

mm3a-mm3c 

2.5 

12000 

mm3b-mm3c 

2.4 

11900 


Uncertainties are ^0.1" or ^490 AU. 


3.2. Far-infrared continuum emission 


Figure |2] presents a comparison of the Herschel PACS 70, 100 
and 160 /urn data with the 3.2 mm continuum data from the PdBI. 
At 70yum, large parts of the central region containing the mil¬ 
limeter cores still appear as a dark extinction silhouette in front 
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Fig. 2. Re-calibrated Herschel PACS images at 70, 100 and 160 n m (from left to right). The scaling in the three panels is from -3.8 
to 4.0, from 5.0 to 17.6, and from 33 to 45 mJy pixel -1 , respectively. The 3.2 mm continuum contours are from 3cr and continuing in 
2<x steps with a ler value of 0.08 mJy beam -1 . A scale-bar and the 3.2 mm synthesized beam are shown in each panel, and the three 
main sources are labeled in the bottom-left panel. The 70, 100, and 160yum PSFs are 5.6", 6.8" and 11.3", respectively. 


of the extended 70yum emission of the hosting molecular cloud 
(Figures [T] & [2]). Although the extinction contrast decreases at 
lOOyurn, we still perceive the central parts as an infrared dark 
cloud (IRDC). At 160/im, the extended emission in the sur¬ 
roundings is not so dominating anymore, and emission structures 
emerge from the inside of the IRDC. 

At the north-western border of the IRDC silhouette, we no¬ 
tice a faint emission point source at 70/im which also persists 
at 100 A/m. In the previou s data reduc t ion pr oducts we used in 
iRagan et all (l2012li and iBeuther et al.l (l2013ah . this source was 
not that apparent. We attribute the differences to the use of 
the newest Scanamorphos algorithms and the inclusion of the 
gyro-correction information (Sect. 12.2b in the current data re¬ 
duction. In order to determine accurate positio ns, we re ferenced 
the PACS maps with the Spitzer/MIPSGAL (iCarev et all 2009 ) 
24 /im data. We did not use individual MIPSGAL maps for refer¬ 
ence, but relied on the recently released MIPS GAL 24 ji m point 
source catalogue ( Gutermuth & Heveti 120151) . Corresponding 
positions in the PACS maps were determined by employing PSF 
photometry using IDL/Starfinder dDiolaiti et alll200dl) . The nec¬ 
essary adjustments were shifts on the order of 1.1-1.6". For 
the previous Herschel data products (with uncorrected point¬ 
ing) we had to use a shift of almost 4" (IBeuther et all l2013al) . 
This supports predictions that the absolute pointing error of 
Herschel data can be brought down to a level of ~079 l<x if 
all the recent correct i ons f or the pointing product are applied 
JSanchez-Portal et al J201 4). This 70 ji m point source is thus lo¬ 
cated very close to mm2. Interestingly, the compact emission 
seen at 160/im is shifted towards the center of the IRDC and 
basically coincides with the mml location (Fig. [2]). The un¬ 
certainties of the peak positions are on the order of 1" for 70 
and 100/tm, due to the position bootstrapping involved. For the 
160 yum astrometry, the uncertainty might even be up to 2", since 
at this wavelength it is hard to find real point sources, and due 
to extended emission and large beam sizes, the effective peak 
position may experience subtle displacements (cf. the compar¬ 
ison of MIPS-70 and PACS-7 0/im p eaks for the EPoS source 
UYSOl reported in iLinz et al.ll2010i) . Still, the formal distance 
between the 70/im point source and the 160/im compact emis¬ 


sion peak is more than 4". Although this is less than the beam 
of the 70/im image, as mentioned in 33.11 source peak positions 
can be determined to much higher accuracy than the nominal 
spatial resolution 0 (down to ^). Hence, we think that there 
is a real shift in peak positions when going to longer wave¬ 
lengths. This may be explainable by the combined action of two 
effects. First, the 160/im data trace colder dust expected in the 
deeper interior of the IRDC, compared to the 70/im wavelength 
range. Second, the 160/im beam is around 1173. It thus con¬ 
volves emission from larger areas, and therefore may comprise 
emission from mml (dominating) and mm2 (minor contribu¬ 
tion), while the finer 70/im beam of 576 can still distinguish 
between mm2 (associated with 70/im emission at its north-west 
side) and mml (without noticeable 70/im emission). 

Table 3. Far-infrared fluxes of mm2 


A 

S 

(jum) 

(mJy) 

70 

90 

100 

357 

160 1 

<339 

3230 2 

0.9 


Flux uncertainties at 70 and 100//m are ~10%. 

1 Upper limit from the residual flux toward mm2 after subtracting the 
point s ource centered on mml. 

2 From IBeuther et al.l ( 2013 a) with an accuracy of ~15%. 

Hence, the PACS data show tentatively that at least object 
mm2 may not be totally starless anymore, but that star forma¬ 
tion processes might have begun at its north-western border. 
Whether the 70 yum point source is already a recently formed pro¬ 
tostar, or just a temperature enhancement created by various pro¬ 
cesses triggered from within mm2, is not easy to answer with the 
presently available data. At 160/im, we eventually see the cold 
dust emission from the bulk of the core material in the mm 1/2 
region. 

The Herschel data can also be used to get a rough estimate 
of the luminosity of the mm2 core. In the 70 and 100yum band, 
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1.07mm continuum N 2 H + [83.0,85.0]km/s N 2 H + [85.0,88.0]km/s N 2 H + [89,92.0]km/s 
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Fig. 3. N2H + (3 - 2) and 1.07 mm continuum data. The left panel and the contours in the three other panels always show the 1.07 mm 
continuum data starting from 3 cr and continuing in 1 cr steps with aim value of 0.6 mJy beam -1 . The color-scale in the 2nd, 3rd and 
4th panels present the N2H + (3 - 2) integrated emission over the velocity regimes marked above each panel. The scaling ranges in 
panel 2 to 4 are [0,0.15], [0,0.25] and [0,0.4] Jy km s -1 , respectively. The synthesized beam of the continuum and line data is shown 
each time in the bottom-left. A scale-bar can be found in the left panel. 


we were able to derive the far-infrared fluxes toward mm2 from 
fitting the point-spread-function (PSF) to the compact emis¬ 
sion source visible in Fig. |2] The uncertainties for the 70 and 
100yum flux measurements are ^10%, including -5% calibra¬ 
tion uncertainty and another -5% from the PSF photometry. 
At 160/mi, this is more difficult because the emission peaks at 
mml. Therefore, in this band we can only derive an upper limit 
for mm2 by fitting a point source to mml, and afterwards de¬ 
riving the mm2 upper limit from the residual image. The 3 mm 
flux density was taken from [Beuther et al. 1 2013a ). Fitting a 
spectral energy distribution to these four data points assuming 
a modified blackbody function accounting for the wavelength- 
dependent emissivity of the dust, we get an estimate of the lumi¬ 
nosity L and cold dust temperature 7d us t. For mm2, this results 
in L - 16^ 4 L© and T dust - 15^ K. The error budget in L and 
r dust includes the uncertainties of the flux calibration as well as 
the selected dust model. While the dust temperature reflects the 
overall cold nature of this region, the comparably low internal 
luminosity of mm2 around 16L© shows that this region with a 
large mass reservoir of -800 M© (section [3711) still has formed 
no high-mass star yet. It should also be noted that the low lu¬ 
minosity cannot be explained well with accretion processes on 
compact protostars because even then the accretion luminosity 
is expected to be higher (e.g., Krum holz et al.ll2007l) . Hence, the 
observed luminosity may likely stem from accretion processes 
on larger surfaces, re-enforcing the early evolutionary stage at 
the onset of star formation with active collapse processes. 

Looking closer at mm3, there is no emission at 70 and 
100//m, and the 160/mi image exhibits only a very weak ex¬ 
tension from mml in the direction of mm3. This is in stark con¬ 
trast to the 3.2 and 1.07 mm data which exhibit almost the same 
fluxes for mml and mm3. This difference is most likely due to 
even lower temperatures within the mm3 region compared to 
mml. 

3.3. Spectral line emission 

Figure [3] presents several N2H + (3-2) emission maps integrated 
each time over different velocity regimes. The velocity regimes 
are selected based on the integrated N2H + (3-2) spectrum pre¬ 
sented in the bottom panel of Fig. [4] We divide the integra¬ 
tions into a low-velocity part between 83 and 85kms -1 , an 


intermediate-velocity regime between 85 and 88kms -1 and a 
high-velocity part from 89 to 92kms -1 . The immediate result 
of that is that the different velocities trace different parts of the 
dense gas. While the low-velocity component is mainly asso¬ 
ciated with the sub-sources mml and mm3, the intermediate- 
velocity component emits toward all three mm cores, and finally 
the high-velocity component emits only toward mm2. Although 
the N2H + (3-2) line has hyperfine-structure as well, the satellite 
lines are comparably weak, and the velocity structure is mainly 
caused by real structure imposed on the main hyperfine compo¬ 
nent. Nevertheless, we also fit the spectra taking into account 
the full hyperfine structure of the line (see below). Independent 
of that, these integrated emission images already show that this 
high-mass starless clump is far from being a kinematically ho¬ 
mogeneous and potentially calm structure, but in contrast to that, 
we see a kinematically complex region that may be at the verge 
of collapse (see also [Ragan et al.1120151) . 

To investigate the kinematics in more depth, we extracted the 
spectra toward the three main mm cores mml, mm2 and mm3, 
which are shown in Figure® Again, the multiple velocity com¬ 
ponents are apparent. To quantify the peak velocities v pea k and 
full-width half maximum Av, we fitted the full N2H + (3-2) hy¬ 
perfine structure to these lines with multiple components. Note 
that, in contrast to the 1-0 transition, for the 3-2 transition of 
N 2 H + the spectrum is much more strongly dominated by the cen¬ 
tral line component. From the 29 hyperfine structure components 
more than 60% of the relative intensities are within the central 6 
components separated by only -0.06 km s -1 whereas the remain¬ 
ing 23 components share the rest of the emission. The result¬ 
ing fits are shown in Table ® One should keep in mind that the 
selection of the number of fitted velocity components is partly 
ambiguous. For example, for mm2 we fitted in total three com¬ 
ponents but one could also have fitted the high-velocity feature 
in mm2 with three components alone and have a fourth compo¬ 
nent around the intermediate-velocity gas. Therefore, these fits 
do not claim to be the final answer, but they nevertheless ade¬ 
quately represent the current data. While the multiple compo¬ 
nents toward individual peaks are interesting in themselves, also 
the derived line width are important. They range in full-width- 
half-maximum (FWHM) between 0.3 and 1.3 km s -1 where the 
0.3 km s -1 have to be considered as an upper limit because with 
our velocity resolution of 0.2 km s -1 no narrower lines can be 
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resolved. The broad end of the distribution is most likely also 
an upper limit, however this time caused by unresolved underly¬ 
ing spectral multiplicity in the lines. With the thermal line width 
at 15 K of -0.15 km s _1 , these new data now resolve the spectra 
into almost thermal lines, just with multiple components. 

For comparison, Figure[4]also shows a single-dish C 18 0(2- 
1) spectrum at 11" resolution (Ragan et al. priv. comm.) as 
well as the PdBI N2H + (3-2) spectrum integrated over all three 
mm continuum cores. While the C 18 0 spectrum easily traces 
the main intermediate-velocity component, however spectrally 
not well resolved in the sub-components (a Gaussian two- 
component fit to the C 18 0(2-1) spectrum results in Av values of 
-1.7 and ~0.3kms _1 , a 1-component fit has Av - 2.0kms _1 ), 
the high- and low-velocity components are not seen at all in 
the C 18 0(2-1) spectrum. The reason for that is likely at least 
twofold and lies within the lower spatial resolution of the single¬ 
dish data as well as the about two orders of magnitude lower 
critical density of C 18 0(2-l) compared to N2H + (3-2) (~ 10 4 
versus - 10 6 cm -3 ). With the new high-resolution observations 
of this dense gas tracer, we are able to really dissect the densest 
portions of this very young high-mass star-forming region (see 
section POl) . 


Table 4. Hyperfine structure fits toward the mm peaks 



^peak 

(km s -1 ) 

Av 
(km s -1 ) 

mml 

83.5 

0.3 

mml 

84.3 

0.3 

mml 

86.3 

0.3 

mm2 

85.8 

1.3 

mm2 

89.9 

0.7 

mm2 

91.0 

0.3 

mm3 

84.4 

0.6 

mm3 

86.2 

0.9 


4. Discussion 

We now have the ability to track the hierarchical fragmentation 
and global collapse of starless gas clumps capturing smaller and 
smaller scales. 

4.1. Fragmentation 

Following llBeuther et ali ( 2013b ), we can estimate the average 
densities p for the large-scale single-dish data as well as the pre¬ 
vious intermediate-scale 3 mm continuum data, corresponding 
to the bottom-left and bottom-middle panels in Figure \T\ The 
estimated average densities p derived from the two datasets are 
2.5 x 10 5 cm -3 and 1.5 x 10 6 cm -3 , respectively. Assuming fur¬ 
thermore average temperatures of 15 K, iBeuther et al.1 (I2013bl) 
estimated the Jeans length Tj ean s and Jeans mass Mj eans that 
predict the expected fragment properties of the correspond¬ 
ing smaller scales in the framework of this isothermal gravita¬ 
tional Jeans fragmentation picture. The Tj ea ns and Mj eans based 
on the large-scale single-dish data are 10000 AU and 0.37 M©, 
whereas the corresponding Tj eans and Mj ean s derived from the 
intermediate^scale 3 mm PdBI data are 4000 AU and 0.15 M 0 
(iBeuther et alll2013bh . 

If we now compare the predicted length and mass scales with 
the corresponding observed values on the different scales, we 
find correspondences as well as differences. Regarding the Jeans 


length and the observed projected separations, the predicted 
Jeans length from the single-dish data of -10000 AU is slightly 
smaller than the core separation in the PdBI 3.2mm data (Table 
0, but the values are consistent within a factor of a few. Going 
to smaller scales, the predicted Jeans length from the 3.2mm 
data of -4000 AU is roughly consistent with the projected sepa¬ 
rations we find again at the smaller scales of the 1.07 mm contin¬ 
uum data (Table [2]). Hence, from a pure length-scale argument, 
the observations of IRDC 18310-4 would be roughly consistent 
with classical thermal Jeans fragmentation. However, our obser¬ 
vations reveal only projected separations, and the real ones can 
be roughly up to a factor 2 larger. 

Does this picture also hold for the Jeans masses? In fact, 
it does not because the predicted Jeans masses on the different 
scales of 0.37 and 0.15 M 0 , respectively, are significantly lower 
than what is found in the 3^2 mm data (masses for mml to mm3 
between 18 and 36 M^. IBeuther et af;2013b ) as well as in our 
new 1.07 mm observations (masses between 2.2 and 13.7 M 0 , 
Table [I]). Hence, while the length scales are roughly consistent 
between the classical Jeans predictions and the data, the masses 
deviate by up to two orders of magnitude. We note that these ob¬ 
served masses are even only lower limits because large fractions 
of the gas are filtered out. 

Similar discrepancies were recently reported by I Wang et afi 
( 2014 ) in their fragmentation study of parts of the “Snake” 
filament. They also found projected separations in their data 
that correspond reasonably well to the estimated Jeans length, 
whereas the fragment masses exceeded the Jeans masses by sim¬ 
ilar margins as in our data. They argued that the typical thermal 
Jeans mass is calculated with the thermal sound speed depend¬ 
ing on the temperature of the gas. However. I Wang et all ( 2014 ) 
evaluate a turbulent Jeans mass using the turbulent velocity dis¬ 
persion of the gas instead of the velocity dispersion based on 
the thermal sound speed. Since the Jeans mass depends on the 
velocity dispersion to the third power, even velocity dispersion 
increases of a factor of a few allow to shift the turbulent Jeans 
masses by more than an order of ma gnitude in their regime of 
observed masses. Hence, [Wang et ak ( 201 4) argue that turbu¬ 
lent Jeans fragmentation may explain the observed properties. A 
first-order inconsistency could arise in this picture when consid¬ 
ering the lengths scales. In the turbulent Jeans scenario, also the 
turbulent Jeans length has to be adapted according to the turbu¬ 
lent velocity dispersion. Although the velocity dispersion only 
enters the equation for the Jeans length linearly, nevertheless, it 
increases the turbulent Jeans length by a factor of a few. And in 
that picture, the projected separations of their cores would be a 
factor of a few smaller than the predicted turbulent Jeans length. 
However, taking into account potential projection effects, this 
difference may be less severe. 

This turbulent Jeans fragmentation picture becomes ques¬ 
tionable for IRDC 18310-4 if one considers that toward one 
of our most massive sub-condensations mml a with 6.6 M 0 the 
measured line widths only shows upper limits of 0.3kms _1 . 
Hence, no significant turbulent contribution to the line-width is 
observed here. 

A different approach to solve the discrepancies between the 
classical Jeans fragmentation and our data is to investigate the 
initial conditions, in particular the initial density structure. While 
the Jeans fragmentatio n analysis starts with uniform and infi¬ 
nite gas structures (e.g. JStahler & Pall5l2005l) . the derived den¬ 
sity structures of the o bserved gas clumps exhibit much steeper 
density profiles (e.g., Beuther et al. I 20021: iMueller et al.l [20021: 
lHatchell & van der Takll2QQ3l) TGirichidi s et al.l (1201 ll) analyzed 
the different fragmentation properties during star formation in 
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numerical simulations varying the turbulent velocity fields as 
well as the density profiles from uniform to Bonnor-Ebert sphere 
and density power-law profiles p oc r~ p with indices p of 1.5 and 
2. As expected, the steeper the initial density profile is, the less 
fragmentation they find and the more massive the final fragment 
masses are. In the extreme cases, the difference in masses rises 
by up to a factor -20. While they do not explicitly give the frag¬ 
ment separations, the fact that they find less fragments, indicates 
that the fragment separation should also increase. 

To summarizes these effects, most likely, the observed frag¬ 
mentation properties in IRDC 18310-4, as well as other sources 
in the literature, could be explained by turbulence-suppressed 
fragmentation of gas clumps with non-uniform density struc¬ 
tures. Magnetic fields can be another agent in inhibiting frag¬ 
mentation fe.g. JCommercon et all201 lb iPillai et aOl2015l) . 

4.2. Kinematic properties 

As outlined in section EU our new high-spatial and high- 
spectral resolution observations of the dense gas tracer N2H + (3- 
2) resolve the spectral lines from the several cores into multi¬ 
ple spectral components where the individual components are 
only slightly above the thermal line width. Although the data 
are not good enough yet to actually resolve coherent ther- 
mally dominated structure s like those found for example in B5 
dPineda et all 120101 12013) . we are reaching a regime for the in¬ 
dividual sub-cores within a high-mass star-forming region that 
is not that far apart from coherence. 

However, in addition to the narrow line widths found for 
individual components, the omnipresence of multiple compo¬ 
nents toward each core implies additional highly dynamic pro¬ 
cesses within the overall high-mass star-forming region. The 
most straightforward interpretation of these multiple velocity 
components is within the picture of a globally collapsing gas 
clump where the different velocity components trace separate 
individual gas parcels that fall toward the gravitational center of 
the whole gas clump. The spectral line N2H + (l-0) signatures 
of such a global collapse were simulated bv [Smith et al . ( 201 3), 
and they find similar spectra with multiple components along 
individual line of sights as seen in our IRDC 18310-4 data. 

The C 18 0(2-1) data allow us also to derive a rough estimate 
of the virial mass. Following M acLaren et al.1 (Il988h assuming a 
density profile p oc 1/r, a radius of the clump of 0.25 pc (Fig. [I]) 
as well as the 1-component fit Av - 2.0kms _1 to the C 18 0(2- 
1) spectrum, the approximate virial mass is -190M©, about a 
factor 4 lower than the gas mass of 800 M© derived from the 
dust continuum data (Beuth er";t al.Ll2013blPl 

In the framework of a globally collapsing gas clump, one can 
use the observed spectral velocity differences also for a simple 
collapse time estimate. Assuming that the difference between the 
velocity peaks is due to cores sitting at different points in a glob¬ 
ally collapsing region, we get converging velocity differences 
along the lines of sight of 2.8 km s -1 in mml, 5.2 km s -1 in mm2, 
and 1.8 km _1 s in mm3. Ignoring for this estimate the plane of the 
sky velocity component that we do not know, these velocity gra¬ 
dients would bring together fragments 10000AU apart in only 
- 1 x 10 4 to - 2.6 x 10 4 yrs. The collapse of the clump may also 
cause the fragments to merge or increase in density, which could 
help explain the hig h frag ment masses in the previous section. 
For example JSmithet aD (l2009l) found, using synthetic interfer¬ 
ometry observations of simulated gas clumps, that the number of 

1 Density profiles steeper than 1/r would result in even lower virial 
mass estimates. 


fragments decreased and their mean column density increased as 
the clumps collapsed. This could efficiently increase the mass of 
the fragments. 

Combining the multiple components with the individual nar¬ 
row line widths as well as the low virial mass, we interprete the 
kinematic data in this region as most likely caused by a dynam¬ 
ical collapse of a large-scale gas clump that caused multiple ve¬ 
locity components along the line of sight, and where at the same 
time the individual infalling gas structures have very low levels 
of internal turbulence. 


5. Conclusions 

Resolving the mm continuum and N2H + (3-2) emission at sub- 
arcsecond resolution (linear scales down to -2500 AU) toward 
the pristine high-mass starless gas clump IRDC 18310-4, we re¬ 
veal the fragmentation and kinematic properties of the dense gas 
at the onset of massive star formation. Zooming through differ¬ 
ent size scales from single-dish data to intermediate- and high- 
angular resolution PdBI observations, the resolved entities al¬ 
ways fragment hierarchically into smaller sub-structures at the 
higher spatial resolution. While the fragment separations are still 
in approximate agreement with thermal Jeans fragmentation, the 
observed core masses are orders of magnitude larger than es¬ 
timated Jeans masses at the given densities and temperatures. 
Hence, additional processes have to be in place. However, tak¬ 
ing into account non-uniform density structures as well as ini¬ 
tial turbulent gas properties, observed core masses and projected 
separations are consistent with cloud formation models. 

While most sub-cores are (far-)infrared dark even at 70 pm, 
the re-reduced Herschel data reveal weak 70 pm emission toward 
core mm2 with a comparably low luminosity of only -16 F©. 
Since such a low luminosity can be caused neither by an internal 
high-mass star nor by strong accretion onto a typical embedded 
protostar, this re-enforces the youth and early evolutionary stage 
of the region. 

The spectral line data reveal multiple velocity components 
with comparably small width (still above thermal) toward the 
individual sub-cores. Relating these data to cloud collapse sim¬ 
ulations, they are agreeing with globally collapsing gas clumps 
where several gas parcels along the line of sight are revealed 
as individual spectral features. The narrow line widths in the 
regime 0.3 to 1 kms -1 indicate that even during the dynamical 
global collapse, individual sub-parcels of gas can have very low 
internal levels of turbulence. 
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Fig. 4. N2H + (3 - 2) spectra toward the three main peaks mml, 
mm2 and mm3 with the multiple component hyperfine-structure 
fits shown in green. The fit parameters are presented in Table 0] 
The bottom spectrum shows a comparison between a single-dish 
C 18 0(2-1) spectrum extracted toward the main peak in grey¬ 
scale, whereas the red spectrum is an average of the three above 
N2H + (3 - 2) spectra (multiplied by 8). The units in the bottom 
spectra are different for both as marked in the side-label. 


































































